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A systematic truncation of the many-body Hilbert space is implemented to study how electrons 
in a quantum dot attached to conducting leads respond to time-dependent biases. The method, 
which we call the dynamical 1/N approach, is first tested in the most unfavorable case, the case 
of spinless fermions (N = 1). We recover the expected behavior, including transient ringing of the 
current in response to an abrupt change of bias. We then apply the approach to the physical case of 
spinning electrons, N — 2, in the Kondo regime for the case of infinite intradot Coulomb repulsion. 
In agreement with previous calculations based on the non-crossing approximation (NCA), we find 
current oscillations associated with transitions between Kondo resonances situated at the Fermi 
levels of each lead. We show that this behavior persists for a more realistic model of semiconducting 
quantum dots in which the Coulomb repulsion is finite. 

PACS numbers: 73.63.Nm, 73.63.Kv, 71.27. +a 



I. INTRODUCTION 

The behavior of strongly interacting electrons confined to low spatial dimensions and driven out of equilibrium is still 
poorly understood. Recent advances in the construction of small quantum dot devices have opened up the possibility 
of studying in a controlled way the nonequilibrium behavior of strongly correlated electrons. For instance the Kondo 
effect, which was first observed in metals with dilute magnetic impurities^, has now been seen in measurements of the 
conductance through a single-electron transistor 2 (SET), in accord with theoretical predictions 3 . Most experiments 
so far have focused on steady state transport through a quantum dot. 

Several theoretical approaches have been developed to calculate electrical transport properties in a far from equi- 
librium situation. Currents generated by a large static bias applied to the leads of a quantum dot in the Kondo 
regime have been analyzed in some detail using the Keldysh formalism combined with the non-crossing approxima- 
tion (NCA)4. Remarkably, stationary currents may also be found exactly with the use of the Bethe ansatz 5 . Exact 
treatments at special points in parameter space have also been carried oul^. More recently, attention has been paid to 
time-dependent phenomena in quantum many-body systems&i&SiifiiiAiiSii^. By looking at the response of an inter- 
acting dot to a time-dependent potential, Plihal, Langreth, and Nordlander studied the several time scales associated 
with different electronic processes^. The response to a sinusoidal AC potentialS*^*^ in the Kondo regime has been 
studied both in the lowi^ and high-frequency limits^. 

NCA has been the most commonly used approximation to obtain response currents due to an external applied 
pulse. The exact solution at the Toulouse point corroborates many of the NCA results^. However, NCA calculations 
have mostly focused on the Kondo regime at not too low temperatures, because in the mixed-valence regime, or at 
low temperatures, NCA is known to give spurious results^. In practice NCA has usually been limited to studies 
in the U — > oo limit. For a typical SET, the ratio of the Coulomb repulsion in the dot to the dot-lead coupling 
strength ranges over [//Aw 10 to 20. Other numerical techniques, such as the Numerical Renormalization Group, 
$_i ' have also been usecU^ to calculate the conductance through dots. However, these methods seem to be limited to 
the static linear-response regime. Perturbative RG methods^ have been used to analyze nonequilibrium transport 
through dots in the Kondo regime^. Recent progress in extending density matrix renormalization group (DMRG) 
methods to explicitly time-dependent problems has been recently achieved with the Time-dependent Density-Matrix 
Renormalization-Group (TDMRG)2i. 

In this paper we introduce a dynamical 1/N approach, where N counts the number of spin components of the 
electron. Physical electrons correspond to N = 2 (spin-up or spin-down). We will also have occasion to consider 
spinless electrons (N — 1) as this case permits a comparison with known exact results for noninteracting electrons. 
Higher values of TV > 2 are of actual physical interest too, as these can occur when there are additional orbital and 
channel degeneracies. The model we study possesses a global SU(N) spin symmetry. 

Unlike NCA, the dynamical 1/N approach is systematic in the sense that it includes all Feynman diagrams up 
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through a given order in 1/N; not only a certain class of them. As crossing diagrams are neglected in the NCA it is of 
interest to compare the two approaches. The static version of the 1/N expansion is a type of configuration-interaction 
(CI) expansion of the sort familiar to quantum chemists. The dynamical 1/N approach has some similarities to a real- 
time perturbation scheme along the Keldysh contour developed by K6nig2£. Perturbing in the powers of the lead-dot 
coupling generates an increasing number of particle- hole excitations in the leads. In practice a resonant-tunneling 
approximation is implemented to restrict the set of diagrams to be considered. Only single particle-hole excitations in 
the leads were included in the off-diagonal elements of the total density matrix, and the intra-dot Coulomb repulsion 
U was taken to be infinite. 

We determine the response currents through a quantum dot under the influence of both a small step bias, and a 
large pulse. Transient oscillatory phenomena are found. In particular, a type of ringing in the current discovered 
previously for non-interacting electrons is recovered within the 1/N approach upon setting N — 1. We compare the 
behavior of the response current for spinning and spinlcss electrons, and demonstrate that the intradot Coulomb 
interaction changes both the period of the oscillations and the decay rate of the currents. We further demonstrate 
that the period of the oscillations, in the case of spinning electrons, does not depend on the dot level energy when 
the dot level is moved from the Kondo into the mixed valence regime, and thus extend previous results in the Kondo 
regime based on NCA&. 

The paper is organized as follows: In Section^we introduce a generalized time-dependent Newns- Anderson model 
for describing a quantum dot and its systematic solution using a 1 /N expansion of the many-body wavefunction. Also 
we describe how non-zero temperature can be treated in the method and how observables are calculated. In Section 
IIIII we find the current through a quantum dot for a small symmetric bias, deep in the linear response regime. In 
Section ll VI we analyze the response of a quantum dot to a large amplitude bias for spinless and interacting electrons. 
We conclude in Section[V]by discussing the advantages and limitations of the dynamical 1/N approach, as well as its 
future prospects. Details of the equations of motion and the calculation of currents are presented in an Appendix. 

II. THEORETICAL APPROACH 

In this section we discuss the correlated-electron model that we employ to study time-dependent currents in a 
quantum dot. We then present its systematic solution order by order in powers of 1/N. 

A. Generalized Newns- Anderson Model 

The model we use to analyze transport through a quantum dot is defined by the following generalized time-dependent 
Newns- Anderson Hamiltonian: 

H(t) = 5>W(*)A + ei 2 )(i)P 2 ] c^c aa +J2^a(t)cf a c km 

a ka 

+ 7= E{^teW^ + ^(*)A]ct^ + fr. c } 

a; ka 

+ ^ E U aan a (n a - 1) + ^ U ab n a n b . (1) 

a a>b 

Here c^ CT is an operator that creates an electron in the dot level a with spin a and c^ a creates an electron in a 
level k in the a-lead. Throughout the paper we adopt the following notation for the indices: Greek letters a, (3, 
7 and 6 label the left and right leads. Greek letter a labels the spin, and a sum over the N spin components is 
implied whenever there is a pair of raised and lowered indices. Operator n a = c\ a c a< T counts the occupancy of 
level a in the dot. Projection operators Pi and P% can be written in terms of n a and project onto the single- and 
double-occupied dot, respectively. Parameters e« (i), ^a\t), V^(t) and V^(t) are respectively the orbital energies 
and hybridization matrix elements for the quantum dot with one or two electrons. The projection operators enable 
the use of different couplings and orbital energies depending on the number of electrons in the dot. For simplicity, 
however, here we assume that the hybridization matrix elements are identical, and independent of level index and 
time: V^} {t) = V^l (t) = V. We also assume that the electronic levels in the dot are the same (apart from the 

charging energy U aa ): eip(t) — £^\t) — e a (t). In Eq. ^ the original hybridization matrix elements are divided 
by a factor of \/N so that the dot level half-width due to the hybridization with either lead, A Q = 77/^1 V^| 2 = A, 
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is independent of N, remaining finite in the N — > oo limit. The density of states per spin channel of either lead is 
denoted by p, here is taken to be constant. Coulomb interaction U aa is the repulsive energy between two electrons 
that occupy the same dot level a, while U a b is the interaction between electrons in different levels. Time-dependence 
may come in through the conducting leads, e ka (t) = e ka + $ka(t), the dot level, e a (t) — e a + &a{t), or both. Functions 
$ a (i) and $ka(t) may have any time dependence. For instance Fig. illustrates the case of a rectangular pulse bias 
potential applied to the left lead keeping both the dot level and right lead unchanged. 

To permit a comparison of the dynamical 1 /N approach to other methods we further simplify the above Hamiltonian. 
In the following we only consider the case of a single level, a = 0, so for notational simplicity we define U aa = U in 
the rest of the paper. The model then corresponds to the usual one considered by many other authors that treats 
only a single s-wave level in the quantum dot. However, we stress that the Hamiltonian, Eq. can also be used to 
study more general and experimentally relevant models. 

B. 1/N expansion of the wavefunction 

The many-body wavefunction is constructed by systematically expanding the Hilbert space into sectors with in- 
creasing number of particle- hole pairs in the leads. Sectors with more and more particle- hole pairs are reduced by 
powers of 1/N in the expansion. The approach was originally introduced to study magnetic impurities in metals2£, 
mixed-valence compounds^ and it was first applied to a dynamical atom-surface scattering problem by Brako and 
Newns^. We have previously applied it to the scattering of alkaliS&Si and alkaline-earth ions such as calcium off metal 
surfaces^ where in the latter case various Kondo effects may be expectedS^Si. The expansion of the time-dependent 
wavefunction for the lead-dot-lead system up to order 0(1 /N 2 ) in the spin-singlet (more generally, SU(N)-singlet) 
sector may be written in terms of the scalar amplitudes f(t), b a . ka (t), e Lljka (t), d ka>q0 (t), s a . Ll ^ a ^ fj {t), a a . >L ^ tkatg0 (t), 
9Li,PS,ka,qf3(t), and hL^,PS,ka,qP\f)- 

!*(*)) = /(*) 1°) + Yj Kka(t) \a;ka) + ^ e LlM {t) |Z/y, ka) + ^ d ka , q p{t) \ka, q/3) 

a; k,a L, fc,7,a g<fc,a,/3 

+ ^2 s a;Lj,ka,qi3(t) \a;Lj,ka,q(3) S + ^ a a - Ll ,ka,qi3(t) \a; Z/y, ka, q/3) A 

a; L, g<fc,7,a,/3 a; L. q<k,~f ,a,/3 

+ ^2 9Lj,P5,ka, q p{t) \Lj,PS,ka,q(3) S + ^ h Ll ,ps,ka, q f3(t) \Lj, PS, ka, qf3) A 

L>P, q<k,y,6,a,0 L>P, q<k,-y,S,a,0 

+ {rest of Hilbert space} . (2) 

The first line of this equation contains the 0(1) amplitudes, the second line has the 0(1/ZV) terms, and the amplitudes 
in the third line are 0(1 /N 2 ). Upper case Roman letters L and P label lead levels above the Fermi level, whereas 
lower case letters k and q label lead levels below the Fermi energy. Expressions for the basis states in the case of 
atoms scattered off metal surfaces can be found elsewhere^, but we rewrite them here for clarity. The vacuum state 
|0) represents an empty dot, with both leads filled with electrons up to the Fermi energy. The remaining states, which 
are all SU(N)-singlets due to the contracted spin indices a and a', are: 

\a;ka) = — L= c\ a c kaa |0) 
V N 

\Lj,ka) = -j= c^c fcoff |0) 



\ka,q(3) = -j==== cl a c kaa cl a c qPry :\Q) 



\a;Lj,ka,q(3) S 



\a;Lj,ka,qP) A 



| Z/y, PS, ka, q(3) 



\L 7 ,P5, ka,q{3) A 



1 



y/2N(N-l) 

1 

y/2N(N + 1) 

1 

y/2N(N-l) 

1 

y/2N(N + 1) 
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We set the Fermi energy to zero (ep = 0) before the pulse is applied (t < U). An ordering convention is imposed on 
the above states such that e q < < and > ep > 0. Fig. [21 shows a schematic representation of the different 
Hilbert space sectors appearing in the l/N expansion up through 0(1/N 2 ). Each row represents a different order 
in the l/N expansion, increasing upon moving downwards in the diagram. The physical interpretation of each of 
the sectors is as follows: State |Z/y, ka) represents an electron excited to level L in lead 7 along with a hole in level 
k in lead a. State \ka, qf3) represents configurations in which two electrons occupy the lowest level (a — 0) of the 
quantum dot simultaneously, with two holes left behind in the leads. These configurations are suppressed in the limit 
U — > 00. States \a;Lj,ka,q(3) S and |a; L7, ka, qf3) A are symmetric (S) and antisymmetric (A) combinations of the 
configuration with one electron in level L and holes in level k of lead a and in level q of lead (3. The division of the 
sector into two parts reflects the fact that the state produced by an electron hopping to the dot from a continuum 
level k while another electron is excited from q to L can be distinguished (because the electrons carry spin) from the 
state in which k and q are interchanged. A similar decomposition is carried out for the sectors with two particle-hole 
pairs described by amplitudes g and h. 

We neglect, at 0(1 /N 2 ), the sector corresponding to a singly-occupied dot with a hole and two particle- hole pairs 
as the amplitude for this sector has 5 continuum indices. This sector is expected to be important for the physically 
interesting case of the Kondo and mixed- valent regimes with near single occupancy of the quantum dot. Its inclusion 
would presumably improve the behavior of the dynamical l/N approach at long times, and we leave this for future 
work. Strictly speaking, through 0(l/iV 2 ) we should also include sectors corresponding to a doubly-occupied quantum 
dot with one or two particle-hole pairs in addition to the two holes in the leads. However, as these configurations are 
described by amplitudes with up to 6 different continuum indices, and as the amplitudes for these sectors is small 
(due to the repulsive Coulomb interaction), the computational work required to include the amplitudes is excessive, 
and we drop them from the equations of motion. 

The equations of motion for the amplitudes that appear in Eq. are given in the Appendix. We note that the 
terms we keep in the l/N expansion is equivalent to summing up Feynman diagrams, including the crossing ones, 
up to order l/N 2 . The inclusion of crossing diagrams is significant because it is these diagrams that are known to 
be responsible for the recovery of Fermi liquid behavior at low temperatures^!, and the disappearance of the Kondo 
effect in the spinless N — > 1 limit. 

C. Calculation of observable quantities 

We conclude our discussion of the dynamical l/N method by explaining how observable quantities such as the 
currents are calculated. At zero temperatures, the initial state of the system prior to application of the bias is chosen 
to be the ground state, which is obtained by the power method. Integrating the equations of motion forward in time 
then yields \^>(t)), and expectation values of observables, O, are calculated periodically during the course of the time 
evolution: 

(6(t)) = <*(t)|(5|*(*)> • (4) 

For the particulars of how to calculate the expectation value of the current operators, see the Appendix. 

Most of the results presented in the paper, with the exception of those shown in the final figure, are for the case of 
zero temperature. At non-zero temperatures we must extend the method^. Suppose that at an initial time t = 0, prior 
to application of the bias, we could find the complete set of energy eigenstates and values {| V I'„(0)), E n } satisfying 
H\^ n (0)) = E n \ty n (0)}. Each of these eigenstates could then be evolved forward in time yielding {(^(i))}- The 
combined thermal and quantum average of any observable quantity, O, would then be given at time t by: 

(OW> = i^(*nW|0|* n (t))exp{ - f^X (5) 

where 

is as usual the partition function of the system. However, as the Hilbert space is enormous in size, the determination 
of the whole set of time-evolved many-body wavefunctions is prohibitively difficult. Instead we sample the Hilbert 
space by creating a finite, random, set of wavefunctions: 71 = 1, ... , n sam pie}- The wavefunctions are generated 

by assigning random numbers to the amplitudes that describe the different sectors of the wavefunction, Eq. [3 Each 
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of these randomly generated wavefunctions is then weighted by half of the usual Boltzmann factor, 

|#„(0))=exp|-^}|H Ti ) , (7) 

and the resulting set of weighted wavefunctions is time-evolved to yield {I'J'n (£))}■ (The wavefunctions are weighted by 
half the usual Boltzmann factor so that thermal expectation values, which depend on the square of the wavefunctions, 
have the correct weight.) Thermal averages of each observable such as the current are calculated periodically during 
the course of the integration forward in time: 

"sample 

{0{t)) = - J2 (*n(*)|0|*„(i)) (8) 

n=l 

where 

Is.mpli Usomple 

Z= £ (* n (i)|*„(t)> = (*n(0)|*n(0)> • (9) 

n— 1 n— 1 

At the very lowest temperatures only a single sample is needed as the approach then reduces to the zero-temperature 
power method described above. At non-zero but moderate temperatures, thermal averages over n samp i e — 30 random 
wavefunctions suffice, and the convergence of the thermal average may be easily checked by increasing the number 
of samples. We stress that time dependence comes in only via the wavefunctions. This makes the present approach 
rather straightforward, avoiding cumbersome Green's function formalisms based on the Keldysh or Kadanoff-Baym 
techniques. On the other hand, the approach is not useful at very high temperatures as an excessive number of 
samples is then required. 

We now present an intuitive picture of how electronic transitions take place as the current flows through the 
quantum dot. Fig. [SJshows some of the sectors, up to 0(1/N 2 ), that contribute to the current through a dot subject 
to the rectangular pulse bias shown in Fig. When the bias is switched on, electronic transitions between the dot 
and the leads take place. Sectors shown in Fig. [21 in which an electron travels from the left to the right leads creating 
cross- lead particle-hole excitations make the most important contribution to the current. The current through the 
dot is a consequence of the formation of particle-hole pairs, with holes accumulating in one lead and particles in the 
other. Sectors of the Hilbert space with increasing numbers of particle-hole pairs, and hence higher order in 1/N, 
become populated as time goes on. Electronic transitions back down to lower order sectors also occur, of course, since 
the Hamiltonian is Hermitian. However, because the phase space of sectors with increasing numbers of particle-hole 
pairs grows rapidly with the number of pairs, when the system is driven out-of-equilibrium, reverse processes back 
down to lower orders occur less frequently. This irreversibility may be quantified in terms of an increasing entropy^. 



III. RESPONSE TO A SMALL, SYMMETRIC, STEP IN THE BIAS POTENTIAL 

We first apply the dynamical 1/N approach to study the response of a quantum dot to a small, symmetric, step 
bias potential. We use parameters appropriate for a semiconducting quantum dot. We either take U — > oo or, more 
realistically, set U — 1 or 2 meV comparable to the values reported by Goldhaber-Gordon et al. for the single electron 
transistor built on the surface of a GaAs/AlGaAs heterostructure that led to the reported Kondo effect 2,32 . We 
further set the dot-lead half-width to be A = 0.4 meV, and e a = — 2A. Thus in the case of spinning electrons the 
dot is in the Kondo regime. In a real SET there are several quantized levels with a typical spacing of order 0.4 meV. 
However, as mentioned in the previous section, here for simplicity we consider only the case of a single energy level in 
the dot. At time t, = 0.5?i/A the bias on the left lead is suddenly turned on, raising the Fermi energy by $/2, while 
the right lead is shifted down by the opposite amount, -$/2, and the dot level is held fixed. The size of the bias is 
chosen to be $ = 0.05 A which is small enough to induce a linear response in the current. 

The leads are assumed to be described by a flat band of constant density of states with all energies taken to be 
relative to the Fermi energy. The band is taken to be symmetrical about the Fermi energy with a half-bandwidth set 
at D = 4 meV. We use M = 30 discrete levels both above, and below, the Fermi energy in our calculations, except for 
the 0(1/N 2 ) amplitudes gLi,P8,ka.,qp{t) and hL~f,ps,ka,q/3(t). As these amplitudes each have four continuum indices, 
we retain only 10 levels above, and below, the Fermi energy. This turns out to suffice as there is little change when 
only 5 levels are retained. To improve the discrete description, the continuum of electronic states in the conducting 
leads is sampled unevenly: the mesh is made finer near the Fermi energy to account for particle-hole excitations of 
low energy. Specifically, the discrete energy levels below the Fermi energy are of the following form: 

e k = [e^-WI M - 1], fc = l,2,...,M (10) 
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with a similar equation for levels above the Fermi energy. Thus, the spacing of the energy levels close to the Fermi 
energy is reduced from the evenly spaced energy interval of D/M by a factor of 7/(e 7 — 1) w 1/13 for a typical choice 

(i) 

of the sampling parameter 7 = 4. The matrix elements V^ a are likewise adjusted to compensate for the uneven 
sampling of the conducting levels 2 ^. We note that the discretization of the continuum of electronic states in the leads 
means that we can only study pulses of duration less than a cutoff time: t cuto ff as 2ttT%/8 : where 8 m D/(13M) is the 
level spacing close to the Fermi energy of the conducting leads. As we show below, however, the 1/N expansion itself 
imposes a more severe restriction on the reliability of the method at long times. 



The case of spinless electrons deserves special attention as it provides a stringent test of the 1/N expansion, and 
is also exactly solvable. For an Anderson impurity model at equilibrium it is known 2 * that the occupancy of the 
impurity at N = 1 is accurate to within 1% of the exact result when terms in the wavefunction expansion are kept 
up through order 1/N 2 . In the dynamical case shown in Fig. |3 however, we find that for e a = — 2A the currents 
decay in time despite the fact that the bias remains turned on. Insight into the breakdown of the dynamical 1/N 
expansion can be gleaned from the equilibrium problem as discussed by Gunnarssson and Schonhammcr 24 . They 
observe that the parameter range e a <C — A is an unfavorable one for the 1/N expansion. The 0(1) approximation in 
this limit differs qualitatively from the 0(1/N) and 0(l/iV 2 ) approximations. In particular there is a large change in 
variational ground state energy (97% in the case of e Q = —A) going from the 0(1) to 0(1/N) approximations. This 
compares to a change of only 18% in energy for the e a = case. In the dynamical problem the observation likewise 
suggests that the case e a = — 2A depicted in Fig. 131 is also unfavorable for the 1/N expansion as higher order sectors 
may be expected to contribute substantially to the wavefunction. By contrast the other N = 1 case of e a = shown 
in Fig. |21 exhibits a clear plateau in the currents. We interpret this behavior as follows: As the dot energy level is 
raised from e a — — 2A to 0, the many-body wavefunction moves closer to the beginning "root" state of the Hilbert 
space expansion, the state |0) in Fig. [21 for which the dot level is empty. Higher-order terms in the 1/N expansion 
are less important in this limit, and the expansion is more accurate. 

The limited number of particle- hole pairs retained in the 1 /N expansion means that electron transfer from one lead 
to the other ceases when the higher order sectors become significantly populated. Therefore, we conclude that the 
dynamical 1/N approach is limited to rather short times after the step in the bias is applied. Inclusion of higher-order 
sectors in the 1/N expansion would presumably permit accurate time-evolution for longer times. The amplitude of the 
0(1/N 2 ) sector for the case of Fig. is quite small when the current starts to decay, consistent with the observation 
that the amplitude of the empty dot sectors (first column in Fig. [2} are small in comparison to the singly-occupied 
sectors (second column of Fig. 0), in the the Kondo and mixed- valence regimes. Including the 0(1 /N 2 ) sector 
consisting of two particle- holes in the leads, an electron in the dot and an extra hole in the leads should improve the 
behavior of the current at longer times. However as a practical matter the rapid growth in computational complexity 
at high orders in the 1/N expansion prohibits extensions to arbitrarily high order. Cluster expansions based upon 
exponentials of particle-hole creation operators 33 may offer a way to obtain more satisfactory behavior at long times. 



Physical electrons have spin, and since N = 2 is a more favorable case from the standpoint of the 1 /N expansion, we 
expect (and find) improved behavior. The Coulomb repulsion U reduces the size of higher-order corrections as it acts 
in part to suppress charge-transfer through the dot that leads to the formation of the particle-hole pairs. Furthermore 
the Kondo effect in the limit e a <C —A is recovered already at 0(1) in the 1/N expansion. As shown in Fig. |21 for 
U — > oo there is a strong enhancement in the current in comparison to the spinless case (N — 1). This is as expected 
from the increase of spectral weight at the Fermi level of the leads due to the Kondo resonance. The current is further 
enhanced at finite U — 1 meV. This can be understood from 1/N calculations of the spectral density at equilibrium. 
A decrease in U results in an increase in the spectral weight at the Fermi energy, reflecting a corresponding increase 
in the Kondo energy scale (see Fig. 9 of Ref. Y2M . The point is also illustrated by an analytical estimate of the Kondo 
temperature for finite U based on weak-coupling renormalization-group scaling24: 



Here Adot — A_l + A# = 2 A is the half- width of the dot level due to its hybridization with both leads. This formula 
yields Tjf w 4K. As the weak-coupling expression Eq. ^2 is only reliable at small values of U, it likely overestimates 



A. Noninteracting spinless electrons (N — 1) 



B. Interacting spinning electrons (N = 2) 




(11) 
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the Kondo temperature in the present case of U/A = 2.5. It is worth comparing the Kondo temperature so obtained 
with the value extracted from the U — > oo Bethe-ansatz expression*, 

k B T K = D — [ — ) ex P(^"-;)' ( 12 ) 



2tt V ttD 



where C = 0.577216 is Euler's constant. Taking N = 2, Eq. [T21 gives only ps 185 mK. 

Hence, the increase in the Kondo temperature obtained from Eqs. ^| and ^2 as t7 is decreased from infinite to 
finite values is consistent with the enhancement of the current shown in Fig. |3| We use Eq. ^]to estimate the Kondo 
temperature at finite-CT in the following section. 



IV. RESPONSE TO A LARGE RECTANGULAR PULSE BIAS POTENTIAL 



In this section we study the response of the quantum lead-dot-lead system to a large rectangular pulse bias potential, 
driving the system well beyond the linear response regime. As the response of the dot-lead system is highly non-linear, 
this situation is quite instructive. We again study both the TV = 1 and TV = 2 cases, and use the same values for e a 
and U as in the previous section but with A = 0.2 meV. The rectangular pulse bias potential is applied as follows. 
At time U = 0.5h/A, a sudden upward shift of the bias with amplitude $ is applied to the left lead, leaving the dot 
energy level unchanged, and the right lead unbiased as depicted in Fig. Q] At a later time t = tf the bias is turned 
off. 



A. Noninteracting spinless electrons (TV = 1) 

The abrupt rise in the bias generates a ringing in the response currents as a result of coherent electronic transitions 
between electrons at the Fermi energy in the leads and the dot level. The period of these oscillations was predicted 
to be±*: 

2tt^ , s 

*• = W^y (13) 

Fig. 0] shows, for the case e a — — 2A, the response current between the left lead and the dot upon setting N = 1 
in the l/N expansion. The oscillation period accords with Eq. ^| For the sake of comparison, Fig. 0|also shows 
the current for the case of interacting, spinning, electrons. Note that oscillation period changes when interactions are 
turned on. As discussed in the next subsection, this is a consequence of the formation of Kondo resonances at the 
Fermi energy in each lead. It is remarkable that the l/N approach reproduces, in the N = 1 limit, the period of the 
oscillations expected from Eq. ^| Unlike NCA, the l/N approach recovers the expected features of non-interacting 
electrons in the N = 1 limit. In particular, many-body Kondo resonances disappear— at N = 1 as they must. As 
mentioned above, the recovery of independent-particle physics can be attributed to the presence of crossing diagrams 
in the l/N expansion. 

Also of interest is the duration of the initial response peak. The transient response time is shorter in the case 
of noninteracting spinless electrons. The transient response fades away more slowly because the repulsive Coulomb 
interaction inhibits electron motion through the dot. 



B. Interacting spinning electrons (TV = 2) 

In equilibrium it is well known that a Kondo resonance can form at the Fermi energy of the leads. For the steady 
state case of constant bias potential, previous NCA calculations of spectral density predicted the splitting of the 
Kondo peak* into two resonances, one at the Fermi energy of each lead. The situation studied here is different as 
we analyze the response to a short pulse rather than the steady state behavior. Nevertheless, the behavior is still 
consistent with the split Kondo peak picture. 

Fig. [3] shows that the response current does not obey Eq. ED Instead the current oscillates with a period that is 
independent of e a , and depends only on the bias. The period can be seen to be 2itTi/<&, which agrees with Eq. ^| 
only upon setting e a — 0. Fig. EJalso demonstrates that the period of the oscillations remains nearly constant as the 
energy of the quantum dot level is varied from e a = (mixed- valence regime) to e a = — 3 A (Kondo regime). Large 
biases split the Kondo peak in two, and the resulting many-body resonances are separated by energy $. Electronic 
transitions between these two peaks induce oscillations of period 27r/i/$. The l/N calculation agrees qualitatively 
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with the NCA calculation of Plihal, Langreth, and Nordlander— . Figs. [S] and [H] show that the split-peak interpretation 
holds even for very large biases, <& = 2A to 10A, and also in the mixed-valence regime. Fig. [5] also shows that the 
magnitude of the response current gradually decreases, and the transient oscillations damp out more quickly, as the 
dot energy level is moved downwards in energy. Insight into this behavior may again be attained from study of a 
quantum dot in equilibrium. In the limit e a / A — ► — oo the Kondo scale vanishes and the spectral weight at the Fermi 
energy is suppressed. The dot enters the Coulomb blockade regime, inhibiting the flow of current. 

Finally we discuss the effect of non-zero temperature on the response currents. Fig. 0shows, for both noninteracting 
spinless and interacting spinning electrons, the currents at three temperatures: T — 100 mK, T — 300 mK and 500 
mK. The latter temperature is sufficiently high that multiple particle-hole pairs will be excited, as the level spacing is 
only D/(13M) as 120 mK at the Fermi energy; therefore the Hilbert space restriction to at most two particle-hole pairs 
is a severe limitation. Nevertheless, there is little change in the response of the spinless electrons as the temperature 
is increased. This is as it should be since the temperatures are still well below the Fermi temperature of the leads. 
By contrast the N = 2 case shows a significant decrease in the current magnitude, and especially the oscillations, 
at the higher temperatures. Using the finite-C/ formula Eq. El with U/A = 10, yields Tk ~ 600 mK. Considering 
the uncertainities involved in estimating the Kondo temperature (Eq. 1111 gives the order of magnitude but the exact 
multiplicative factor is unknown) we can expect a suppression of the Kondo resonance and its associated effects at 
500 mK. The behavior of the currents shown in Fig. 0is consistent with this interpretation, as the period of the 
transient oscillations decreases to approach that of the spinless system. 

V. CONCLUSIONS 

We have presented a method for calculating the response of a quantum dot to time-dependent bias potentials. 
The approach, which is based upon a truncation of the Hilbert space, is systematic because corrections can be 
incorporated by including higher powers in the 1/N expansion. In agreement with previous approaches we find 
coherent oscillations in the response currents at low temperatures. We note that although the frequency of these 
oscillations is of order a terahertz for parameters typical of quantum dot devices, it should be possible to detect the 
oscillations experimentally—. The dynamical 1/N approach permits a realistic description of the quantum dot as 
finite Coulomb repulsion and multiple dot levels may be modeled. For typical device parameters, response currents 
are qualitatively similar to those found in the U — - > oo limit. However, the magnitude of the current is larger. 

We also discussed an extension of the dynamical 1/N expansion to treat the case of non-zero temperature. In 
the case of interacting spinning electrons the magnitude of the currents, and the nature of the transient response 
are sensitive to relative size of the system temperature in comparison to the Kondo temperature. In contrast non- 
interacting spinless electrons show little temperature dependence. 

The dynamical 1/N approach complements other methods such as NCA and the TDMRG algorithm. The primary 
limitation of the dynamical 1/N approach comes from the truncation of the Hilbert space in which only a small 
number of particle-hole pairs (in this paper, at most two) are permitted. Hence, as it stands the method is not well 
suited for the study of steady state situations, and we cannot use it, for instance, to answer the question of whether 
or not a lead-dot-lead system exhibits coherent oscillations at long times^Si. As NCA sums up an infinite set of 
Feynman diagrams, including ones with arbitrary numbers of particle-hole pairs, it provides a far better description of 
long-time behavior, including steady state current flow. For the same reason NCA is also superior at high temperature. 
However finite Coulomb interactions and other realistic features of lead-dot-lcad systems are technically difficult to 
incorporate within NCA, but not in the dynamical 1/N method. Also, NCA shows unphysical behavior in the mixed 
valence regime, at very low temperatures, and in the spinless N = 1 limit. The dynamical 1/N approach does not 
suffer from these pathologies. The method shares some features with the TDMRG algorithm as both are real-time 
approaches that truncate the Hilbert space in a systematic, though in a different, fashion. Unlike the present approach, 
TDMRG can treat strong electron interactions between electrons in the leads and can be used to study tunneling 
between Luttinger liquids^. On the other hand, non-interacting electrons in the leads are treated exactly within the 
dynamical 1/N approach. In TDMRG these require as much or more computational effort as interacting electrons. 

A complete description of the electronic transport properties in a SET would require taking into account the whole 
energy spectrum of the quantum dot and would include both the direct Coulomb repulsion and spin-exchange between 
electrons in different energy levels of the dot. Also the conducting leads should be described by realistic densities 
of states. The dynamical 1/N approach is well suited to accommodate these complications. Another aspect worth 
further attention is the case of orbital degeneracy of the dot levels. As the degeneracy increases, the Kondo resonance 
strengthens, and its width decreases^. A combined experimental and theoretical exploration of the effect of orbital 
degeneracy on transport properties would be interesting. Finally, it should also be possible to extend the method to 
describe two coupled quantum dots attached to leadsSi. Competition between Kondo resonances and spin-exchange 
between the two dots leads to unusual features such as a non- Fermi liquid fixed poinfeiSiii. Probing nonlinear transport 



9 



as the parameters are tuned through this fixed point 42 could possibly shed some light on the physics of heavy-fermions 
materials. 
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VI. APPENDIX 

In this Appendix we present details of the equations of motion and the calculation of the currents. 

A. Equations of Motion 

Equations of motion appropriate for the atom-surface scattering problem have been published elsewhere^ 7 -. We 
rewrite them here, generalizing the equations to include two sectors that are of order 0(l/iV 2 ) and substituting two 
conducting leads for the metallic surface. To remove diagonal terms in the equations of motion we introduce the 
phase factor 

0«(t) = i f Q ef{t')dt' , (14) 
which is the phase of the quantum dot level when it is decoupled from the leads, and 



1 r* 

ka(t) = ^ J e ka (t')dt' , (15) 



the phase of the electronic levels in the decoupled leads. Upon projecting the time-dependent Schrodinger equation 
onto the different sectors of the Hilbert space, the equations of motion may then be extracted. Following Refs. |2g and 
12^ we use capital letters to denote amplitudes in which the diagonal phases of the corresponding lower case amplitudes 
have been factored out: 

ih J t F = E V £& exp{# fca (*) - <f> { a\t)}} B a , ka 

a;ka 

ihj t B a , ka = V$ a exp{#«(*)-<fe a (i)]}F 



+ S afi ^1-1/NJ2 V© exp{-i[t/ - <M*) + 2 ^o 2 V) - $\t)]} M k - ?) D k^ qf 3 + % - k) D qP . ka ] 
+ "7= E V £lr ex P{#i 1} W ~ <M*)]} E LlM 



a;q/3 



\l E V aZ e Mil<t> q p(t) - 4> ( aHt)]W(k - q) A a , Lj ^ qf3 - 6{q - k) A a . <L ^ M ] 

a;qj3 



ih-D ka , qP = y/l - l/N exp{i[U - cj> qP (t) + 24> ( z \t) - (t)]} B Q . ko 
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+ y/l - l/N V$ a exp{i\U ~ <M«) + 2^\t) - $\t)]} B 0;qfi 

+ \IJjY, V & exp{i[U-<f> L7 (t)+2$ ) (t)-$\t)}} S QMa , qP 
hi 

= 6 a , \f^V$; exp{-i{U-cp Ll + 2$ ) (t)-$\t)}}D ka ^ 

+ \J^T exp ^ i ^ 1) W ~ ^ m E ^ ka + V SL «P{i[^ 1} (*) ~ M)}} E Lj , qP ] 

+ ~m S exp{i[<Mi) ~ ^ (*)]}[*(£ ~ J)G a ;In,J6MAf> + d ( J ~ L)G a , J& ^, ka ^] 



ih^A a . MM>q p = exp{i[0«(i) - 9/9 (t)]} £ L7 , fctt - V$ a exp{#«(t) - te (i)]} 

+ E exp{i[0j,(t) - - ■/;//.,:/ -. /,>„,,< - 0(J - L)^ o; j 5ji7 , fc a, g /3] 

^ JS 

i^G LjiJStka ,qi3 = -y=v£}g exp{i[4:js(t) - <j>£< ] (t)]}S a ., Ly , ka ,qp + --Lv^* exp{i[0 L7 (t) - <j>W (t)]}S a -j 5 , ka , q 

i^ t H Ll ,JS,k a , q = -y==v£)g enp{i[ct>jg(t) - ffl (t)]}A a , Ll ^ ql3 - -^=V^exp{i[<j) Ll (t) - ^ (t)]}A a . J5M , q p ■ 

Here 9(k — q) and similar terms are unit step functions that enforces the ordering convention on the continuum indices. 
The equations of motion are numerically integrated forward in time using a fourth-order Runge-Kutta algorithm with 
adaptive time steps. We monitor the normalization of the wavefunction to ensure that any departure from unitary 
evolution remains small, less than 10~ 6 . 

B. Calculation of Currents 

A consideration of the rate of change of the quantum dot occupancy, as determined from the commutator [n a , H], 
permits the straightforward identification of the current operators between the dot and the left and right leads: 

a;k ** ^ 

= -^=Im{ 5>$;A + V$:P 2 ]ctc aa } (17) 

where as before the Greek index a refers to either the left or right lead. The expectation value of the current is 
calculated as the quantum and thermal average of the operator Eq. 1171 For each time-evolved wavefunctions, \P n (i), 
the corresponding expected current may be expressed in terms of the amplitudes in the different sectors: 

(*„(t)|Ja|*„(t)) = ™ Imj J2 V $a exp{i[<p ka (t) - <f>M(t)]}F*B a . M 

^ a;k 

+ y/l- l/N Y, VqZ eM~*[Ut - ha(t) + 2^\t) - $\t)]}B*. qf3 [6(k - q)D ka , q p + 6(q - k)D q0>ka ] 



k,q/3 



a;L,qf3 



+ 9(q — k)S a - : Lj, qa ,kl3) 



2N 

a;q,L-y.kp 



a;q,Z/7,fc/3 

+ S F S exp{-*[(7t - <M*) + 2$\t) - ^\t)]}S* Q , Lakl ^D kl , qP 

L,kf,q/3 
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4>a (t)]}Sa-,jp,krf,qs[0(L - J) G la,jf3,k~/,qS + ®( J ~ L ) G Jf3,La,k~/,qs} 

Note that the current remains finite in the N — > oo limit as the amplitudes appearing in Eq. 1181 are at most 0(1). 
Thus Eq. 1181 is the total current summed over all N spin channels and divided by a factor of JV, and hence equivalent 
to the current per spin channel (since spin-rotational invariance remains unbroken) . Fig. [5] depicts some of the 
contributions to the current between the dot and the leads. At each time step we check that current is conserved 
to within an accuracy of 10 -6 (when charging of the dot is taken into account). A thermal average of the current 
operator, using Eq. [H] is the final step: 

(J Q (*)> = |£(#„(i)|Ja|*„(i)>- (19) 



+ 7f E v^Mi[<PMt) 

* a;L,J/3,kf,q5 
V a;L,J3,k-y,qS 
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FIG. 1: Energy level representation of interacting electrons in a quantum dot before, during, and after application of a bias 
pulse of amplitude $ to the left lead, holding both the right lead chemical potential and the dot energy level, e a , fixed. Here 
U is the Coulomb interaction between electrons inside the dot. 
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FIG. 2: Schematic representation of Hilbert space as organized by the 1/N expansion. The diagram shows the succesive 
electronic transitions due both to the coupling of the dot levels to the leads and the applied bias. Particle-hole excitations are 
produced in both leads. Sectors with increasing number of particle-hole pairs appear at successive orders in the 1/N expansion. 
We retain configurations which are at most of order 0(1/N 2 ). Not shown, but included in the many-body wavefunction, arc 
configurations with two particle-hole pairs in a single lead, or two particles in one lead and two holes in the other lead. 
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FIG. 3: Response of the electric current per spin channel between the left lead and a quantum dot to a small symmetric bias 
step potential. Time is given in units of h/A where A = 0.4 meV, and current is in units of e/h. A small, symmetric, step bias 
potential of amplitude $/2 = 0.025A is applied to the left lead and an opposing bias of —0.025 A is applied to the right lead. 
The arrow at time t = U = 0.5H/A indicates the moment when the step in the bias is turned on. The electronic level in the 
dot is held fixed at energy e a = — 2A. (For comparison in the N = 1 case we also show results for e a = 0.) The non-interacting 
spinless case as calculated within the dynamical 1/N approximation upon setting N = 1 is compared to the case of interacting 
electrons (N = 2) with U = 1 and 2 meV and U = oo. Currents initially grow and then decay (evident in the time range shown 
here only for the case N = 1 and e a = — 2A) despite the fact that the bias remains on. The eventual decay of the currents 
is a consequence of continued production of particle-hole pairs, saturating the high-order sectors of the 1/N expansion. The 
response current is strongly enhanced in the interacting case with respect to the spinless case due to the presence of a Kondo 
resonance at the Fermi level of the leads. For (7 = 1 and 2 meV, the current is further enhanced reflecting an increase in the 
Kondo temperature. 
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FIG. 4: Electric current per spin channel between the left lead and a quantum dot due to application of a large rectangular 
pulse bias potential. Time is given in units of ft/A where A = 0.2 meV, and current is in units of e/h. Here the bias of 
amplitude $ = 5A is applied only to the left lead (see Fig. 0, and the level in the dot is held fixed at e a — — 2A. The bias is 
turned on abruptly at time t = 0.5H/A and off at time t = 3.5H/A. In the same plot we compare the non-interacting, spinless 
case (N = 1) with the N = 2 case for both the U — > oo limit and for the experimentally relevant value of U = 2 meV of a 
SET—. The transient response time of the interacting electrons is longer than for non-interacting electrons. As can be seen 
upon comparison with the time scales (horizontal lines) the period of the oscillations also increases as N changes from 1 to 2, 
in accord with the formula t p = i DU t only if we set e a = in the case of spinning electrons. 




FIG. 5: Dependence of the electronic current per spin channel between the left lead and the quantum dot on the amplitude, 
$, of the externally applied rectangular pulse bias potential (see Fig. 0. Time is in units of h/A where A = 0.2 meV, and 
current is in units of e/h. We consider interacting spinning electrons (N = 2) and take U — > oo. The quantum dot level is 
either in the Kondo regime (e a /A = —2) or in the mixed valent regime (e a /A = —0.5). The pulse starts at time t = 0.5R/A 
and ends at time t — 3.5h/A. The current is plotted for different applied biases: $ = 2A, 5A, and 10A. A comparison of time 
scales (horizontal lines) shows that the period of the oscillations does not obey Eq. 1131 but rather t p = 2nh/<&. 
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FIG. 6: Electric current per spin channel between the left lead and the quantum dot for different dot level energies. Time is 
in units of ft/A where A = 0.2 meV, and current is in units of e/ft. We lixJV = 2 and U = 2 meV, values appropriate for 
a semiconducting quantum dot. A rectangular pulse bias potential of amplitude $ = 5A is applied with the dot level set at: 
e a = — 3A, — 2A, and — 0.5A. The pulse is abruptly turned on at time t = 0.5ft/A and off at time t — 3.5ft/ A as marked by the 
arrows. We find that the period of the oscillations remains constant even as the quantum dot level is moved from the Kondo 
into the mixed- valence regime. As e a /A becomes more negative the Kondo scale decreases and the magnitude of the current 
is also suppressed. 
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FIG. 7: Temperature dependence of the response current per spin channel through a quantum dot. Noninteracting spinless 
(N = 1) electrons are compared to interacting spinful electrons (N = 2). Time is in units of h/A where A = 0.2 meV, and 
current is in units of e/h. A rectangular pulse bias potential of amplitude $ = 5A is applied. We set e a = — 2A, and for the 
interacting case we set U = 2 meV. In contrast to the noninteracting spinless case, for interacting spinful electrons the period 
of the oscillations and the magnitude of the current decrease significantly as the temperature is increased, approaching those 
of the spinless system. 



